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ABSTRACT 

The Fermi Large Area Telescope has recently discovered two giant gamma-ray 
bubbles which extend north and south of the Galactic center with diameters and 
heights of the order of if ~ 10 kpc. We suggest that the periodic star capture 
processes by the Galactic supermassive black hole Sgr A*, with a capture rate of 
''"cap ~ 3 X 10"^ yr"^ and an energy release of ~ 3 x 10^^ erg per capture, can 
result in hot plasma injecting into the Galactic halo at a wind velocity of u ~ 10^ 
cm s~^. The periodic injection of hot plasma can produce a series of shocks. 
Energetic protons in the bubble are re-accelerated when they interact with these 
shocks. We show that for energy larger than E > 10^^ eV, the acceleration process 
can be better described by the stochastic second-order Fermi acceleration. 

We propose that hadronic cosmic rays (CRs) within the "knee" of the observed 
CR spectrum are produced by Galactic supernova remnants distributed in the 
Galactic disk. Re-acceleration of these particles in the Fermi Bubble produces 
CRs beyond the knee. With a mean CR diffusion coefficient in this energy range 
in the bubble Db ~ 3 x 10'^° cm^ s~^, we can reproduce the spectral index of 
the spectrum beyond the knee and within. The conversion efficiency from shock 
energy of the bubble into CR energy is about 10%. This model provides a natural 
explanation of the observed CR flux, spectral indices, and matching of spectra 
at the knee. 

Subject headings: acceleration of particles - galaxies: jets - Galaxy: - shock waves 
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INTRODUCTION 



The recent discovery of a pair of giant Fermi Bubbles in the Galactic center (GC) is 
one of the most remark able events in astro physics. The first indication of the structure 
appeared in a paper by 



Dobler et al. 



( 120101 ). which they called the Fermi Haze. Using the 



special p r ocedu re of background subtraction of the Fermi Large Area Telescope (LAT) data, 



Su et al. 



( I2OIOI ) discovered a pair of symmetric structures above and below the Galactic 
plane in the GC direction. The origin of the bubble, if its existence can be proved, is still 
enigmatic, and up to now a few models have been presented in the literature. The team 
that subtracted this structured gamma-ray emission from the total diffuse Galactic emission 
presented different explanations of the phenomenon, but they seemed to favor the model 
that a single massive release of energy in the GC when a huge cloud of gas or a star cluster 
was captured by the central supermassive black hole about 10 Myr ago. 



A similar explanation was suggested by 



Guo fc Mathewsl (120111 ) and 



Guo et al. 



(1201 ih . 



They assumed that the Fermi Bubbles were created by a recent active galactic nuclei jet 
activity about 1-2 Myr ago, which was active for a duration of ~ 0.1-0.5 Myr, releasing 
energy totaling ~ (1-8) x 10^^ erg. The bipolar jets were ejected into the Galactic halo along 
the symmetric axis perpendicular to the Galactic plane. 

It is important to note that the existence of the bubbles wa s first evidenced in X-rays 



detec ted by ROSAT as a narrow envelope with very sharp edges (IBland-Hawthorn fc Cohen 



2OO3I ) and later the Wilkinson Microwave Anisotropy Probe ( WM AP) detected an 



excess of radio signa . 



Dobler fc Finkbeiner 



s at t 



2008 



le location o f the gamma-ray bubbles ([Finkbeiner 



Dobler 



2004; 



2OIII ). The i^O^^T structure is explained clS db fast 
wind with a velocity u^, ~ 10^ cm s~^ driving a shock into the halo gas. This phenomenon 
requires an energy release of about 10^^ erg at the GC and this activity should be periodic 



on a timescale of the order of 10 Myr. This requirement of energy release in the GC is 
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consistent with the observations of the existence of hot plasma with a temperature of about 
10 keV in the GC region but with a radius of 30-50 pc onl y. This cannot be co nfined and 



will escape from the GC with a speed of ~ 10® cm s ^ (IKoyama et al. 



20071 ). Therefore, 

sources with a power of about lO'*^ erg s"^ mus t have heated the plas r na or released ~ 10^^ 
erg from the GC in the past 10 Myr. However, ICrocker fc AharonianI ( 1201 ll ) have proposed 



a relatively slow energy release (~ 10^^ erg s~^) from supernova (SN) explosions as a source 
of proton production in the GC. The observed gamma rays come from hadronic processes 
of the protons in the halo. The plasma in the halo is extremely turbulent and the protons 
are trapped for a time comparable to the Hubble time. But this model requires a separate 
origin of electrons, which have a much shorter lifetime than protons, to explain the WMAP 
data. 

Intensive energy release has been observed, indeed, at the center of normal galaxies 
as strong variations of X-ray radiation. There are common characteristics of these X-ray 
sources. First, all of them have been bright sources, and their X-ray luminosity could go up 
to about 10^^ erg s~^. Second, they have shown a high level of variability in their X-ray 
light curve within years. In the "high state" , the luminosity of one source could be at least 
100 times higher than the luminosity in its "low state" . Third, most of them have a super 



soft spectrum during 



;he flare, with effective bla ckbody temperatures of only about 10-100 



eV f lKomossa fc Bade 



1999 



Halpern et al 



20041 ). The classic examples that satisfy these 



characteristics are RX J1624.9+7554, RX J1242.6-1119A, RX J1420+5334, RX J1331-3243, 
and NGC 5905. Many scenarios have been proposed to explain these phenomena, but 

most of them fail to explain some of the observed results. A detailed discussion of these 

I 1 

scenarios can be found in iKomossa fc Bade! (119991 ). Among all of the listed models, the 



tidal disruption model is the most commonly accepted, and it gives the most satisfactory 
explanation to the observations by considering the radiation from the disk. In this model, 
when a star passes by a black hole within a capture radius, where the black hole tidal 
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force becomes stronger than the self-gravity of the star, the star can be captured. The 



detailed capt ure and disruption process o : 



authors (e.g., 



Rees 



1988 



Cannizzo et al 



a m ain-sequence star has been studied by several 



19901). The capture rate of main -sequence stars 



in our Galaxy and in other galaxies is about 10 ^ yr ^ to 10 ^ yr ^ (see 



1999 




Alexander '. 


2005 


). Recent 


y more stel 


Esquej et al. 


2008 


: Gezari et al. 


2008, 


2009: 



Syer fc Ulmer 



ar capture events have been observed (e 



Komossa fc Bade 



1999 



Cappelluti et al. 



20091). 



Dynamical studies of nearby galaxies suggest that most, if not all, galaxies with a bulge 



component host a cen tral sup ermassive black hole and that ; 



are tightly correlated ( iMagorrian et al 



1998 



Tremaine et al. 



le bu 



2002 



ge and black hole m ass 



Greene fc Ho 



20071 ). 



The Burst Alert Telescope on boa rd Swift has identified a transient X-ray source calle d 



GRB110328A 



Leloudas et al. 



Cummings et al. 



20 111 ) with later optical identification (jCenko et al. 



2011 



201 ll ). which is located in the direction of the constellation Draco when 



it erupted in a series of X-ray flares. The distance to this source is determined to be 



z ~ 0.35 by using and Oiii emission lines by the Gemini telescope ( iLevan et al 



2011b[ ). 



The c haracteristics of GRB 110328A appear inconsistent with those of a gamma-ray 



burst (IBarthelmy et al. 



201ll ). In fact, its time-dependent characteristics including various 



timescales in light curve, multi- wavelengths, etc., seem to b e better explained in terms o f 



the \ 


:;i( 


lal disruption of a star by a supermassive black hole ( 


A 


meida 


2011: 


Bloom et al. 


2011 


Levan et al. 


2011a: 


Burrows et. al. 


2011; 


Zauderer et al. 


2011) 


. All these recent 



observations suggest that stellar capture processes are quite common in other normal 
galaxies. 

Observations have also revealed much evidence of unusual processes occurring 
in the central region of our Galaxy, for instance, the enigmatic 511 keV annihilation 
emission discovered by the International Gamma-Ray Astrophysics Laboratory (see, e.g.. 



Knoedlseder et al 



20051 ) whose origin is still debated. The hot plasma has a temperature of 
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about 10 keV which cannot be confined in the GC and , therefore, sources w ith a power of 
about 10^^ erg s~^ are required to heat the plasma (see 



Koyama et al 



2007. and references 



therein). In fact, plasma outfiow s with velocities of > 10'' cm s ^ are observed from the 



nuclear region of our 



(IBogdan fc Gilfanov 



Galaxy (see 



Crocker et al. 



20101 ) and from the nucleus of Andromeda 



20101) • Time variations of the 6.4 keV line and X-ray continuum 



emission observed in the direction of molecular clouds in the GC which ar e supposec 



a refiection of a g;iant X-ray fiare oc curred several hundred years ago (see 



Ponti et al. 


2010; 


Terrier et al. 


2010 



Inui et al. 



to be 



2009 



and references therein). HESS observations of the GC 



in the TeV energy range indicated an explosive injection of cosmic rays (CRs) there, whic h 



might be associated with the supermassive black hole Sgr A* (e.g., 



Aharonian et al. 



20061). 



In a series of papers ( ICheng et al. 



2006 



2007 



Dogiel et al. 



2009a 



20111), we 



developed a model of energy release in the GC due to star accretion onto the central black 
hole for the interpretation of X-ray and gamma-ray emission from the GC. Our goal was to 
explain these observational data in the framework of a single model. Basic assumptions in 
these models are (1) the Galactic supermassive black hole Sgr A* can capture a star at a 
rate of Us ~ 10~^-10~^ yr"^, and (2) the energy release from each capture in the form of 
a fiux of subrelativistic protons is ~ 4 x 10^^ M^R~^Ml[^{b /O.l)'"^ erg, where M* (in 
units of Mq) and i?* (in units of Rq) are the mass and the radius of the captured star, 
Mbh (in units of IO^Mq) is the mass of the supermassive black h ole, and b is the ratio of 



Alexander 



20051). In a 



the periapse distance to the tidal radius Rt (see the review of 
time scale much longer than the capture timescale, this model can be treated to have an 
average power injection W ^ 3 x 10^° erg s~^. These protons heat the surrounding plasma 
by Coulomb losses to 10 KeV. 



Based on this model. 



Cheng et al. 



(2011 



CCDKI model) argued that up to several 



hundred capture events might have occured in the past 10 Myr, which may have generated 
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a series of shocks propagating through the central part of the Galactic halo and thus 
produced accelerated relativistic electrons responsible for the bubble emission. Processes 
of charged particle acceleration by the bubble shocks in terms of sizes of the envelope, 
maximum energy of accelerated particles, etc., may differ significantly from those obtained 
for SNe. In this paper, we examine whether a "signal" from charged particles accelerated 
in the bubble region can be seen in the spectrum of CRs observed at the Earth. We 
present simple estimations of hadronic CR acceleration by the Fermi Bubble shocks up to 
energies above 10^^ eV. The paper is organized as follows. In Section |2l we review current 
understanding of CR acceleration by supernova remnants (SNRs) and conclude that this 
process can only produce CRs with energies less than 10^^ eV. We present a simple solution 
of the multiple-shock structure in the halo in Section [3l In Section |4] we discuss the protons 
accelerated by the bubble shocks. We emphasize that a broken power law of particle 
distribution must be formed because of the finite spacing between consecutive shocks and 
the spectral break naturally occurs at 10^^(u/10^cms~^)(/sh/30pc)(i?/5;uG) eV. Charged 
particles below and above this critical energy are accelerated by two different acceleration 
mechanisms. In Section El we calculate the total particle spectrum by summing up the 
contribution from all shocks in the bubble and compare it with the observed hadronic CR 
spectrum with energies larger than 10^^ eV. In Section [6l we suggest a model that can 
produce the CR spectrum within and beyond the "knee" (around 3 x 10^^ eV). Summary 
and discussion is presented in Section [3 



2. CR ACCELERATION BY SNRs IN THE GALAXY 



From a general point of view, SN explosions are enough to supply the power needed 



for the luminosity o. 



Berezhko et al 



1994 



CRs in our Ga laxy, Lqr ~ 10"^^ erg s ^ (for a g eneral review se e 



Reynolds 



20081 ). Diffusive shock acceleration (see 



Krymskii 



1977 



Bell 
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19781 ) is considered to be a viable and natural mechanism for CR accelerated by SNRs. 
The mechanism produces a power-law spectrum with the necessary spectral index that is 
observed experimentally. The simplest one-dimensional kinetic equation describing this 
process has the form 



df ^ d fr^\r \ ^ du{z) d . 3j.\ _g 

dt dz \ ^ ^ Qz ) 3p2 dz dp^^ ' 



where z is the coordinate perpendicular to the shock front, p is the particle momentum, 
u{z) is the velocity distribution which describes a velocity jump at the shock, and D is 
the spatial diffusion coefficient. The solution of this equation is a power-law function, 
fip) oc p~"' 1 in which the spectral index 7 is a function of the velocity jump at the shock. 
For strong shocks with a Mach number much larger than unity, 7 = 4. The corresponding 
energy spectral index is z/ = 2 {N{E) oc E~^). 

The current status of the observations of middle-aged SNRs by LAT on broad 



Fermi with an energy range from .2 to 100 GeV 



shock-acceleration theory of SNRs (ICastro fc Slane 



las provided sorne insi ght into the 



2010 



Uchivama 



2011 



]. Assuming 



Castro fc Slanel (120 lOf ) deduced the 



the gamma rays are produced by hadronic processes, 
spectral index of four SNRs ranging between 2.1 and 2.4. However, whether the observed 
GeV gamma rays are produced by hadronic processes or leptonic processes is very difficult 
to differentiate. On the other hand the ambiguity can be removed if broadband emissions 
are observed. In particular if GeV and TeV spectra can be described by a single power law, 
which is steeper than -E"^, the hadronic processes could be the more favorable mechanisms. 
Currently about 10 SNRs have been detected in both GeV and TeV bands, including 



Tycho and CTB37A, whose GeV- TeV gamma-ray emissio n shows 



Caprioli 



unifor mly steep spectral 



20111 ). All these recent 



indices of about 2.3 and 2.2, respectively (see Table 1 of 
observations are consistent with conventional SN shock-acceleration theories, which suggest 
that the spectrum of CRs is roughly described by E~'^. 
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Many fundamental questions related to the assumption that SNRs are the sources of 
Galactic CRs are still open. The maximum energy of the accelerated particles is the main 
concern for this scenario, which can be roughly estimated from a very simple relation. The 
acceleration time at the shock is racc(-E') ~ D{E)/u'^^, where Ush is the shock velocity (~ a 
few 10^ cm s~^). The minimum value of the diffusion coefficient at the shock follows from 
the Bohm diffusion scenario, i.e., -DBohm(-E') = {c/3)rL{E) where is the Larmor radius 
of the particle. Equating the acceleration time with the lifetime of the shock T we get an 
estimate for the maximum energy of accelerated particles after time T, 



-E'max ~ 



Ze(3shUshBT . 



(2) 



where /5sh = "Ush/c, B is the magnetic field strength at the shock. The combination UghB/c 
can be interpreted as an effective electric field. 

For an SNR of typical age tsnr ~ 1000 yr, the maximum energy of protons is 
easily estimated by requi ring t hat the acceleration tim e remain smaller than tsnr- 



Lagage fc Cesarskyi (119831 ) and iBerezhko fc Volkl (120001 ) demonstrated that the maximum 



energy of protons within the scenario of Bo 



for standard galactic SNRs. 



i m diff usion is as large as E 

max 

lO^^-lO^^ eV 



Berezhko et al 



(1l994l ) estimated the efficiency of acceleration 



when a feedback reaction of accelerated particles on the front structure was included and 
they showed that in the Bohm limit CRs absorb about 20% of the explosion energy. The 
acceleration process acts as an effective viscosity in widening the region of the shock velocity 
jump and eventually the acceleration process stops. 

However, outside the quasi-linear model the acceleration of CRs at t h e sho ck fronts 



of S NRs may make th e acceleration of particles more effective. As 



also 



Bykov et al. 



Belli (l2004h (see 



20091 ) showed that during acceleration at shocks of SNR, the magnetic 



non-resonant fluctuations were strongly driven. A nonlinear MHD simulation indicated 
that CR-excited turbulence could amplify the magnetic fleld. It appears that acceleration 
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to the spectral break at 10^^ eV normally ceases as an SNR enters the Sedov phase. Thus, 
CR acceleration by SN shocks can only provide particles with energies less than 10^^ eV. 

The spectral index of the observed CR flux changes from 2.7 to 3.1 around energy 10^^ 
eV, and this is known as the knee. The standard model of CR acceleration by SN shocks 
cannot explain CR energies above the knee, because it only produces a single power-law 
spectrum up to the energy around 10^^ eV. In addition, the CR spectrum flattens again for 
energies above 10^*^ eV, and this is known as the ankle. Large size is required to accelerate 
and to confine charged particles above the ankle (the Larmor radius at these energies is 
comparable with the halo height). The origin of CRs above the ankle is generally attributed 
to an extragalactic origin because those particles could not be confined inside the Galaxy 
and known potential galactic accelerators could hardly accelerate particles to such high 
energies. 

The origin of the steepening for E > 10^^ eV is still an open question and different 
mechanisms o f CR a cceleration in the range 10^^-10^^ eV in the Galaxy have been proposed. 



Ptuskin et al 



(|2010[ ) assumed that the CR flux with energies above 10^^ eV is produced 
by very young galactic SNRs. They modeled the particle acceleration by spherical shocks 
with back-reaction of CR pressure on the shock structure. The signiflcant magnetic fleld 
ampliflcation in young SNRs produced by CR streaming instability may lead to a flux of 
CRs with the maximum energy of accelerated particles about 5 x 10^^ eV. In this model, 
the steepening of the CR spectrum at the knee position is due to distortion of the spectrum 
ejected from young SNRs by the propagation process. 



Anothe r interpretation was suggested by 



torn , and 



Lagutin et al. 



Erlvkin fc Wolfendalel toO& l 



Erlykin et al. 



(120081 ) who a. s sumed that CRs at the knee were produced by 



a single, recent local SN. Recently 



Butt 



(120091 ) summarized problems (including energies 



within the knee) of the conception that isolated SNRs are the main accelerators of CRs and 
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discussed alternative scenarios of CR acceleration. 



The steepening of the CR spectrum at high energies may also be the result of a 



change in properti es of diffusion in t 



mentioned first by 



le interstellar medium. This effect of propagation was 



Syrovatskiil (119711 ) who noticed that the standard diffusion of CRs in 



the interstellar medium might be changed by convection due to a drif t of these part i cles in 



Ptuskin et al 



f ll993h 



the large-scale Galactic magnetic fields. This model was developed in 
who assumed two types of CR diffusion in the Galaxy: the usual diffusion due to particle 
scattering on fluctuations of random magnetic fields and the Hall anisotropic diffusion (drift 
motion) due to the large-scale Galactic magnetic field whose effect might become important 
just above the knee energy. 



As an alternative model, 



Jokipii fc Morfilll (119871 ) suggested a mechanism of acceleration 



in the Galaxy of ultra high energy CRs in a Galactic wind and its hypothetical termination 
shock. In this scenario, SNRs accelerate the bulk of CRs up to 10^^ eV. These particles 
are further accelerated up to 10^^-10 '^° eV at a te r minat ion shock which is at a distance 



of a few hundred kpc from the disk. 



Ip fc AxfordI (I1992I ) analyzed multiple interactions of 



particles with SNRs in the Galactic disk as a source for CR acceleration above the knee. 
liowever, too many shocks are required in the disk in order to produce CR flux at the knee. 



Bykov fc Toptygin 



(I1993I ) showed that regions of CR acceleration to energies above 10 eV 



might be OB associations where concentration of shock fronts is very high. We will discuss 
this model in Section HI 

In summary, it is generally agreed that SNR shocks can only accelerate particles to 
energies less than 10^^ eV. On the other hand, accretion processes in the GC may generate 
giant shocks which are effective for particle acceleration above 10^^ eV. It is interesting to 
note that the notion of recurrent activities in the Galactic center as a source of cosmic 



rays has been proposed some 30 years ago (IKhazan fc Ptuskin 



1977 



Said et al 



1981 
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Ptuskin fc Khazarull98ll ). although not in the same perspective as the present paper. 



STRUCTURE OF SHOCKS IN THE FERMI BUBBLE 



As was assumed in CCDKI the central massive black hole captures a star every 
Tcap ~ 3 X 10^ yr; as a result, about ~ 3 x 10^^ erg of energy in the form of subrelativistic 
particles is released. This heats up the surrounding gas in the central region of our Galaxy. 
The hot gas expands into the halo and forms a propagating upward shock. The situation 



medium (see, e.g.. 



Weaver et al. 


1977; 


Bisnovatvi-Koean & Silich 


1995 


)• 





The gas distribution in the disk and in the halo was derived in 
and can be presented as a double exponential distribution 

2 



Cordes et al. 



fll99lh 



ra(p, z) 



0.025 exp 
+ 0.2exp 



1 kpc 

z 

0.15 kpc 



X exp 



P 



20 kpc 



X exp 



p — 4 kpc 
2 kpc 



cm 



(3) 



The energy release in the GC as a result of star capture can be either impulsive or 
continuous depending on the characteristic times of star capture and energy dissipation of 
subrelativistic protons (plasma heating by Coulomb losses). The capture time is roughly 
Tcap ~ 3 X 10^ yr and the dissipation time Tdiss in the CCDKI model is determined by the 
rate of ionization losses of protons injected with energy which is given by 

10^ 



Tdiss 



n 



En 



.Icm-V V 100 MeV 
where n is t he gas density in t he vicinity of the GC which can be quite high (see the 



(4) 



discussion in 



Cheng et al 



20071 ). e.g., in nearby molecular clouds n > 10^ cm 



If Tcap Tdiss we have the case of stationary energy injection in the GC. In this 
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case the region of heated gas is bounded by a single shock (see I Weaver et aLlll977( ). For 
Tcap ^ Tdiss a muhi-shock structure is formed in the halo with shocks of different ages. 
A similar multi-shock structure can also be created if there are epochs of high-frequency 
star captures in the GC. Thus, the number of shocks is determined by the injection and 
dissipation parameters. 

For a highly simplified case of the exponential atmosphere with the scale height zq, i.e., 
the plasma density n(z) = nn ex pf— z/zn), an analytic s olution of shock propagation wa s 



Bisnovatyi-Kogan fc Siliclj|l995[ ) 



obtained by iKompaneetsI (119601 ) (see also the review of 
This solution gives a qualitative picture of the shock propagation and parameters of the 
medium bounded by the shock that roughly described the situation expected in the Fermi 
Bubble. If the rate of energy injection is L then the radius of the shock as a function of the 
height z and the time t is 



p{z, t) = 2zq arccos <( ^ exp ^ j 



1 



exp 



z 



(5) 



Here z is the coordinate perpendicular to the Galactic plane. 



y 



XLt 



V{t)mno 



1/2 



dt. 



(6) 



V{t) is the current volume enveloped by the shock 

a{t) 

V{t) = 27r y p'^{z,t) dz. 



a{t) is the position of the top of the shock 

a{t) = -2zoln ( 1 



(7) 



2zo 

L = W/Tcap is the average luminosity of the central source, ■jp is the polytropic coefficient, 
and A describes the fraction of explosion energy converted into the thermal energy of gas 
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Fig. 1. — Multi-shock structure in the bubble from the Kompaneets' solution showing seven 
representative shocks. 
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(see iBisnovatyi-Kogan fc Silichlll995[ ). As it follows from Equations (|5]) and ([8]), for the 
finite time ti determined from the condition y{ti) = 2zq the shock breaks through the 
exponential atmosphere and the bubble top a(ti) tends to infinity while the bubble radius 
in the Galactic plane {z = 0) tends asymptotically to p = 2^0 cos~^(l/2) ^ and for 
z zq to p hzq. 

A numerical solution of the system dS])-® in dimensionless coordinates for shock waves 
at different ages is shown in Figure [H This figure is meant to be illustrative only. In reality 
the distribution of shocks should be far more complicated. Neglecting the thin 4 kpc torus 
component of Equation the scale height of the "atmosphere" is = 1 kpc. The shock 
distribution in Figured] is suggestive of a cylindrical bubble with an edge at a radius ps — 3 
kpc. 

Based on this, we put forward a quasi- stationary model of the Fermi Bubble which 
we regard as the source of hadronic CRs with energy larger than 10^^ eV. The essence 
of the CCDKI model is that energy is quasi-periodically injected into the halo when the 
stellar capture processes take place and may exist over a timescale comparable with the age 
of Milky Way. Consequently, the Fermi Bubble should have a stationary structure. The 
idealized Kompaneets' solution above shows that there is a stationary sideway boundary for 
shocks. For quasi-periodic star capture, the bubble interior is filled with shocks propagating 
in series and eventually stopping at pb — 3 kpc. However, in the Kompaneets' solution 
dissipation processes are ignored. The realistic situation has to be described by a set 
of dissipative hydrodynamic equations, which takes account the shock propagation in a 
non-uniform medium and various dissipation processes, including shock heating, energy 
transfer into CRs, slowing dow n due to accumulat ion of material, etc. Fitting with the 



observed gamma-ray spectrum. 



Cheng et al. 



(120111 ) concluded that electrons should have 



an escape time scale of 15 Myr in order to explain the spectral break position. The 
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characteristic dissipation time scale of the shocks should be of the same order if these 
electrons are transported away by the shocks. The shocks will be mostly dissipated when 
they arrive at the sideway boundary (ps) and there will be no pile-up of shocks at ps- As 
the speed of shocks along the bubble axis is progressively larger than the sideway speed, 
the upper and lower boundaries of the bubble will be the same as the halo boundary 
{z = ±H = ±10 kpc from the mid-plane). Thus, the bubble has a stationary structure. 

The same result for the dimension of th e bubble can be obtained if we use the swept-up 
mass model proposed by 



Cheng et al. 



( 120111 ). In this model the radius of the sideway shock 



front or the swept-up front is given by ps = ^2\W/7imnAzu'^p, where Ugp is the speed of 
the sideway shock front (or the swept-up front) and Az = UTc^p- With u ~ 10® cm s""*^, 
Tcap ~ 3 X 10^ yr, Az ~ 30 pc. We argue that the sideway shock front or the swept-up front 
will disappear when its speed is smaller than the local sound speed. In Cheng et al. (2011) 
we have estimated that the shocks heat up the halo to ~ 1 keV and the characteristic sound 
speed is of the order of ~ 3 x 10'' cm s~^. Thus putting Ugp ~ w^, the sideway boundary 
Pb ^ 3.2kpc(AVr/2 x 10^"^ eTgy/'^{Az/30pcy^^^{n/10^^ cm^^y^/'^{vs/3 x lO^cms-i)^^ 
(note that XW is the fraction of injected energy converted into thermal energy of gas). 



4. PROTON ACCELERATION BY THE BUBBLE SHOCKS 

Correct analysis of shock acceleration in the bubble requires sophisticated calculations 
in each stage of this process which we perform latter. Now we present simple estimates of 
the characteristics of the spectra of the accelerated particle in the framework of the CCDKI 
model. 

Below we analyze the spectrum of protons accelerated in the bubble and discuss 
whether the bubble's contribution to the total flux of CRs in the Galaxy may explain the 
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knee steepening. We remind the reader that the generally accepted point of view is that the 
flux of relatively low energy CRs (< 10^^ eV) is generated by SNRs, which eject a power-law 
spectrum E~'^ into the interstellar medium. This spectrum is steepened by propagation 



(escape) pr ocesses in the Ga. 



details, see 



Berezinskii et al. 



axy i n accordance with the spectrum observed near Earth (for 



1990 



chap. 5, Section 6). However these sources can hardly 
produce CRs with energies > 10^^ eV, at which a steepening (the knee) in the CR spectrum 
is observed. For character istics of the knee spectrum and models of its origin see the review 



of iKotera fc Olintol fl201ll ). We suggest that the bubble could generate the flux of CRs at 
energies > 10^^ eV because the shocks in the bubble have much larger length scales and 
longer lifetimes in comparison with those in SNRs. 

In the framework of CCDKI, the bubble may flll with hundreds of shocks propagating 
in series one after another, though a single shock structure cannot be excluded. The average 
separation between two shocks is given by 



^sh — Tc&p'U — 30 



' cap 



U 



pc. 



(9) 



3 X lO^yry VlOScms-i 

However, the exact amount of time between two consecutive shocks depends on the actual 
time between two consecutive capture events and their energy releases. There is another 
important spatial scale which characterizes processes of particle acceleration by a single 
shock: the diffusion length scale at a single shock lo ~ D/u. Here u is the shock velocity 
and D is the spatial diffusion coefficient of the energetic particles near a shock which 
depends on particle interaction with small-scale magnetic fluctuations. In the Bohm limit, 
D ~ crL{E)/3, where Tl^E) = E/ZeB is the particle Larmor radius. In this case 

cr^ cE 



D 



u ZeBu 



(10) 



The problem of particle acceleration in conditions of supersonic turbulence (mul tiple- 



shock structure) has been extensively analyzed (e.g.. 



Spruit 



1988 



Achterberg 



19901 : 
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Schneider 



1993 



Melrose fc Pope 



2003). In a series o 



Bykov fc Toptygin 



pape rs by 



19931) as well as quasi-periodic flows (e.g., 



Bvkov fc Toptvginl fll993f ) 



Webb et al. 



Bvkov fc Fleishmanl fll992[ ) and 



(120011 ) the idea was applied to acceleration processes in OB associations, 
which is quite similar to the structure of the bubble. They introduced a dimensionless 
parameter characterizing the acceleration regimes 



ul 



lb = ~ 

^ Id D 



sh 



ul. 



sh 



The critical energy Ei that separate two regimes of acceleration can be estimated from the 
condition ^/^ ~ 1 or l£,{Ei) ~ /sh- For the conditions of the Fermi Bubble the critical energy 
is 



El 



ZeBul 



sh 



10^^ z 



B 



sh 



U 



108 



cm 



eV. 



(12) 



Bvkov fc Toptvginl fll993[ ) and 



5 /iG / v30 pc^ 

In the case of 1 or ^ ^h, the analysis in 

Bykov fc Fleishmanl fjl992l ) showed that there is a combined effect of a fast particle 
acceleration by a single shock, which generates the spectrum E^"^ and relatively slow 
transformation of this spectrum due to interaction with other shocks (stochastic Fermi 
acceleration) into a hard E~^ spectrum in the intershock medium at relatively low energies. 
However, it is unclear if such slow transformation can be completed within the lifetime of 
the shocks in the bubble. A detailed numerical analysis is needed. Furthermore, from the 
general point of view the characteristic acceleration time is quite short in the range E < Ei, 
which is roughly given by the shock acceleration time ctl/u"^. In the range E > Ei, the 
acceleration time scale increases to the time of stochastic acceleration, c/sh/w^- With an 
average Galactic spatial diffusion coefficient Dq outside the bubble and a Galactic halo of 
height H, the characteristic escape time is Tcsc ~ H^/Dq- We expect that escape processes, 
which play a crucial role in determining the particle spectrum shown in our next analysis, 
are insignificant in the range E < Ei. Therefore the particle spectrum produced by the 
bubble should be ~ E~'^ for E < 10^^ eV, where 2 > u > 1. As discussed in Section [2l 
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SNRs are the major contributors for CRs with energies E < 10^^ eV, the exact particle 
spectrum generated from the bubble is unimportant in the energy range of E < 10^^ eV. 

In the case of ^ 1 or ^ /sh, inside the region of supersonic turbulence the 



acceleration regime shifts to a p ure stochastic acce^ 



extend the equation derived by 



eration by the supersonic turbulence. We 



Bykov fc Toptygiru (119931 ) for CR acceleration in supersonic 



turbulence in a stationary state and axisymmetric geometry to include spatial dependent 
diffusion coefficient and external source, 

where p and z are the cylindrical spatial coordinates and p is the particle momentum. 
D{p,p) is the spatial diffusion coefficient and K,{p,p) is the momentum diffusion coefficient. 
Their spatial and momentum dependence in our model is described below. Q{p, z,p) is the 
possible CR source which will be useful in our numerical example in Section |6l 

As we mentioned above, it is reasonable to assume that CRs with energies E < 10^^ 
eV are supplied by SNRs. Therefore, in this section, we concentrate on the analysis of the 
acceleration of CRs in the energy range E > Ei in the bubble by supersonic turbulence. In 
Section |6l we will treat the case with SNRs and bubble and deal with energy from less than 
10^2 eV to larger than 10^^ eV. 

We set the boundary condition of the distribution function at the Galactic halo outer 
boundary 

/|s = 0, at p = pG and z = ±H . (14) 

Proton acceleration in the bubble depends sensitively on the acceleration parameters and 
structure of the bubble. In the following we present a detailed analysis. We model the 
bubble region as a cylinder extending above and below the Galactic plane from z = to 
z = ±H with a radius p = Pb and assuming there is no CR source inside the bubble (i.e.. 
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(5 = 0). The diffusion coefficients inside and outside the bubble are supposed to be different 



D{p) = Db d{pB -p) + Dg 9{p - pb) 
K,{p,p) = kbp'^0{pb - p) , 



(15) 
(16) 



where Db ^ c/sh/3 is the coefficient inside the bubble as a result of interactions with a 
su personic turbulen c e and Dq is the average diffusion coefficient in the Galaxy, e.g., defined 



m 



Berezinskii et al. 



(119901 chap. 3, Section 6). The momentum diffusion coefficient is 



Kb ~ /Db- The momentum dependence of / is represented by a power-law function, 
f{p) P~"' 1 where 7 should be determined from Equation f|T3|) . 

To understand the dependence of 7 on other parameters, we make two simplifications 

ue of 7 significantly. First, for H < pc, 



of Equation (fT3|l , which do not affect the va. 

' 3 



as expected from 



boundaries at p = pc are insignificant (see 



Strong fc Moskalenkd (119981). particles that esc ape through the radial 



Berezinskii et al 



1990 



, chap. 3, Section 3), and 



we can shift the halo boundary to infinity, i.e., pc = 00. Second, we model the axisymmetric 
geometry of the problem as planar geometry (i.e., we assume df /dp ^ f/p)- We go back 
to the axisymmetric geometry afterward. 



As in iBulanov et al.l (119721 ) and iBulanov fc Dogell (119741 ) we search for solutions to 
Equation (IT^ by the method of separation variables, / = R{p)Z{z)p~^ . The solution for 
Z{z) has a very simple form 



Zniz) 



cos{knZ / H) 



(17) 



where /c„ = 7r(n + 1/2). We should point out that to make sure / is non- negative we must 
take n = for physical solutions. 

To illustrate ideas, we consider the case Dq = Db and approximate the axisymmetric 
geometry as planar (i.e., cPR/dp^ ^ 1/ pdR/dp). Using the dimensionless variable q = p/H, 
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Equation (fT3|) can be simplified as 



^n + 7(3-7) 



Dn 



R = 



(1^ 



which has the exact form of the Schrodinger equation for a rectangular potential well: 

2m 



- f/o) ^ = , for < p < ps 



and 



d'^^ 2m ^ 

-j^ + = , for p> pB 



(19) 



(20) 



when we define £ and f/o as 



£ = — and f/o = — 7 3-7— — 
2m 2m iJ^ 



The solution of this equation is well-known (e.g., 
22): 



Landau fc Lifshitz 



1991 



(21) 

, chap. 3, Section 



^{q) = Ci exp (kng) , for < , where A;„ = V—2mS/h 

"^{g) = C2exp {-kng) , for ^ > , (22) 

^(f?) = sin {qg + S) , for < g < gs , where <^ = \/2'm{E — UQ)/h 



From the continuity of the logarithmic derivative of \i/ at the well boundaries we have the 
condition 

arcsin ( = iil^^lEA (23) 

W-2mUoJ 2 ^ ^ 

and 

cos ^ = ±x^ , for odd j , (24) 

sin{ = ±x^ , for even j , (25) 



where 



^ KpB J h I -2 

^=—— and x = — 1 



(26) 
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For strong acceleration, u^H'^/Dl^ ^ 1, Equation (12^ and f l25|) determine a finite 
number of levels. From these equations, the first term of the series gives 




4 p%nB 



(27) 



m 



For weak acceler a tion, w^ff^Z-Dn ^ 1, we can use the shallow well solution presented 



Landau fc Lifshitzl (jl99ll ). In this case there is only one level at i?o — f^o that gives 



7 




> 1. 



(28) 



A rough estima te of the power of CR pr oduction by the bubble can be done in the same 



way as presented in 



Berezinskii et al. 



(Il990l . chap. 1, Section 4) for Ge V CRs. The energ y 



density of CRs at = 3 x 10^^ eV is uck ^ 6.7 x 10 erg cm 



17 



see 



Kotera fc Olinto 



201ll ). Then the power required for the bubble to produce the knee at the Earth is 

cuc^Mh 



B 



X 



(29) 



where Mh is the total mass of hydrogen in our Galaxy, which is about 10^^ g, and 
X is an extrapolation of energies > lO^feV of the CR grammage derived from the 



chemical composition by 



Jones et al. 



( 1200 ll ) up to energies of about several hundred GeV, 
x{E) ~ 11.8 X (4.9 GeV/-E')°'^^ g cm~^. Then we obtain the required power of CR sources 
at the knee energy range Wb ~ 2 x 10"^^ erg s~^, which can easily be supplied by star 
capture processes. 

More accurate values of 7 can be derived from numerical calculation of the axisymmetric 



case: 



Id ( dRi 
g dg \ dg 



fcn + 7(3-7) 



D 



B 



Ri = 0. 



(30) 



Here the index i = 1,2 denotes the regions inside {g < gB) and outside {g > gB) the bubble, 
respectively. Note that kbi = hb and k,b2 = (see Equation (fT6|) ). A solution inside and 
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10"= 



■Dg/Db = 1 

-Dg/Db = 10 



7 = 1.5+ V9/4 + 7r^DB/KBg^ 



^- 10 




10" 



10 



10" 



10" 



10' 



Fig. 2. — Spectral index for different acceleration rates. 

outside the bubble is searched as series of the Bessel functions {Ju, inside the bubble) and 
the McDonald functions {Ky, outside the bubble). 



The boundary conditions at the bubble radius, Q = Qb are 



Ri{qb) = R2{,Qb) and Db 



dRi 



dg 



D. 



dRo 



G 



Qb 



dg 



(31) 



BB 



These relations can be satisfied if 



ai = kl + 7(3 - ^)kbH''/Db < 0, 



(32) 



and the above requirement implies 7 > 3. For n = we have 



Ri{q) = CiJoiV-ctiQ) and R2{g) = C2Ko{7rg/2) 



(33) 



or 



DbV-(^iMV-(^iQb) iTDGKi{7igB/2) 



2Ko{ngB/2) 



(34) 
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Fig. 3. — Spatial distribution of protons in the central region of the Galaxy as a function 
of Q. Here q is normalized to H (the height of the halo), and in this example the bubble 
boundary is taken as qb = 0.6. 

Numerical results of different ratios of Dq/Db as a function of the ratio between the 
times of particle escape from the Galaxy H^/Db and the time of Fermi acceleration 
I/kb ~ Db/u^ are shown in Figure [2J The dashed curve shows the approximate estimate 
given by Equation ( 128|) . which indicates that Equation (128|) is a good approximation and 
we will use it in the next section of data fitting for simplicity. 

The spatial distribution of accelerated protons for different ratios Dq/Db is shown in 
Figure |3l We can see that these two distributions do not have a qualitative difference even 
if the ratios differ by a factor of 10. 
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5. BUBBLE CONTRIBUTION TO CRs AT ENERGY LARGER THAN 10^^ 

eV 

In this section, we use the general models developed in the last section to fit the CR 
spectrum for E > 10^^ eV but we ignore any possible spectral modulation effect as a result 
of propagation from the bubble to the Earth. We have shown that the charged particles in 
the bubble can be described by a power-law spectrum, 

where Ei is given by Equation ( IT2l) and z/ = 7 — 2, which can be expressed as 



for El < E < E2 where E2 is the high-energy cutoff and will be discussed below. We have 
used the momentum diffusion coefficient kb = v?/ Db- Here A^^o is a normalization constant 
which is chosen to fit the observed CR spectrum by assuming that the CR spectrum 
with energy E > 10^^ eV is entirely contributed by particles from the bubble. E2 can be 
estimated by Equation (jj]), which gives 

{E^Un^SxlO'' (^] (—^] (—^ X eV, (37) 

^ ' VS/iCy VlOMyry V108cms-i/ ' ^ ' 

where T is the time taken by a single shock propagating from the disk to the top of the 
bubble. However, this estimate assumes that particles are only accelerated by a single 
shock. As described in the previous section, the high-energy particles are accelerated by 
multiple shocks in the bubble; in other words, particles can diffuse downward and continue 
to be accelerated by younger shocks in the lower part of the bubble. Therefore T ~ 10 Myr 
is only a minimum lifetime, which means that the above estimate can be considered as the 
lower limit of E2. Another possible way to restrict the maximum value of E2 is when the 
Larmor radius of the particles {tl = E/eB) is larger than the radius of the bubble {H/2), 

10- (^)(^)eV. (38) 
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In fact this estimate should be more appropriate for E2 and we will use it in our model 
fitting process. 

In our model particle spectrum there are four parameters: the lower cutoff Ei 
(Equation (fT2|) ). the upper cutoff E2 (Equation ( 138|) ). the spectral index u (Equation ( 136|) ). 
and the normalization Nq (Equation fl35p ). We consider that Igh = 30 pc, B = 5 fiG, 
u = 10^ cm s~^ and H = 10 kpc are very reasonable mean values in the bubble, therefore in 
our model fitting we fix Ei = 10^^ eV and E2 = 10^^ eV. On the other hand, the conversion 
efficiency from the shock energy into particle energy and the mean diffusion coefficient Db 
in the bubble are the most uncertain parameters. Therefore in our model fitting we treat 
these two parameters as fitting parameters. The solid line in Figure H] indicates the best 
fit. The best fit gives u = 3.12, which corresponds to Db ~ 3 x 10'^° cm^ s~^. This is 
in agreement with the estimation Db ^ c/sh/3 if we take the average separation between 
shocks as 4h ~ 100 pc (see Equation ([9])). This seemingly large diffusion coefficient is in 
fact at le ast an order of m agnitude smaller than the coefficient in the halo for E > 5 x 10^^ 

suggested Dq ~ 2.0 x 10^8(^/4.9 GeV)°-^^ cm^ s'^). 



eV (e.g.. 



Jones et al. 



2001 



If the CRs in the energy range 10^^ eV < E < 10^^ eV arriving at the Earth come 
from the bubble, then the power provided by the bubble for CR in this energy range is 
given by Wcr ~ J^°,''"^ 47iR^EcR{E)dE ~ 3 x 1039(i?/10 kpc)^ erg s-\ where Fcr{E) 
is the observed CR energy flux and R is the mean distance to the bubble. We find that 
the conversion efficient from shock power, W ~ W^/Tcap ~ 3 x 10^*^ erg s~^, is about 10 %, 



which is consistent with recent estimation by using the Fermi-LAT data (lAbdo et al. 
2OIOI ). Realistically, particles could escape from the bubble through various locations of 
the bubble's surface, where the local strength of the magnetic field may be different. In 
addition, we have pointed out that there are two possible ways to estimate the value of 
E2 (Equations fl371) and fl38|) ). Therefore it is likely that Ei and E2 should have some 
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distribution coming out from the bubble. In the dotted hne in Figure S] we keep the spectral 
index as that of the solid line but assume that Ei and E2 have a uniform distribution 
between 2 x 10^'' eV < < 2 x 10^^ eV (which corresponds to 1 yuG < i? < 10 fiG) and 
(-£"2)111111 < E2 < (-E'2)max, respectively. We can see that the model curve starts to drop 
around 3 x 10^^ eV and matches the data better. Furthermore, we have argued that the 
injected spectrum from the shock is E~^ for E < Ei and slowly becomes harder for E~^ due 
to interactions with multiple shocks after a sufficiently long time. However, it is unclear if 
this spectral hardening process can be completed within the finite lifetime of the particles 
in the bubble (see Section H] for a more detailed discussion). Therefore we cannot predict 
the exact spectral index for E < Ei. For reference purposes, in Figure H] we show these two 
possibilities: the dashed line for E'"^ and the dash-dotted line for E~'^, respectively. Both 
El and E2 of these two fines also assume uniform distributions. 



6. MODEL FOR HIGH-ENERGY CRs WITHIN AND BEYOND THE 

KNEE 

The observed high-energy CR spectrum is a broken power law. The energy spectral 
index for CR energy smaller than the knee (around 3 x 10^^ eV) is 2.7. The index increases 
to 3.1 for larger energy (until around 10^^ eV where the spectrum becomes harder again). 
The general argument is that the sources (i.e., acceleration mechanisms and/or sites) 
responsible for the energy range within and beyond the knee should be different. The 
intriguing fact is that the two power laws match quite well at the knee. This coincident 
problem is difficult to solve if the two sources are totally unrelated. 

CRs within (i.e., energy less than) the knee are generally attributed to SNRs in 
our Galaxy (see Section [2]). In the previous two sections, we alluded to the acceleration 
site of CRs beyond the knee (i.e., energy larger than the knee) to the Fermi Bubble. In 
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e co ntribution to the flux of CRs. The data are summarized in 
Kotera fc Ohntol fl201ll ). Th e spectrum has been m ultiphed by E"^ for clarity. 



iments include Tibe t AS-'y f|Amenoniori et al 



KASCAD 
CHiR es-I: 



Grande flApel et al. 



Abbasi et al. 



20081) KASCADE flKampert et al 



yxper 



2004h . 



20091 ). High Resolution F 



20091 ). HiRes-II flAbbasi et al. 



y's Eye cosmi c-ray detector-I 



20081 ). and Auger ( lAbraham et al. 



2010l ). The solid line shows the contribution from the bubble predicted by Equation (!36|) for 
El = 10^^ eV < E < E2 = 10^^ eV. The dotted line has the same spectral index as that of 
the solid line but both Ei and E2 are assumed to obey a uniform distribution (see the text 
for explanation). The dashed line and the dash-dotted line have the same spectral index for 
El < E < E2 and distribution of Ei and E2 as that of the dotted line except the spectral 
index changes to 1 and 2, respectively, for E < Ei (see the text for explanation). 
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Section HI we studied Equation (fT3|) without a source in the bubble. The solution gave the 
characteristic spectrum of the system. In reality we need a source or a seed population for 
the bubble. We deem that we do not need a new seed population for bubble acceleration. 
Instead we propose the following model: some of the CRs produced by SNRs in the Galactic 
disk are re-accelerated in the bubble to energy beyond the knee. This model has the 
potential to solve the coincident problem naturally, because the source of CRs beyond the 
knee is seeded by the source within the knee. 

To consolidate our idea, we work out a concrete numerical model. Essentially, we solve 
the stationary state CR transport equation f[T51) in our Galaxy with two Fermi Bubbles 
(one on each side of the Galactic plane). We modeled our Galactic halo as a cylinder of 
radius pc = 20 kpc, and the top and bottom at ±10 kpc from the mid-plane. Each Fermi 
Bubble is also a cylinder of the same height ±10 kpc, but with a radius Pb = 3 kpc. 

The Fermi Bubbles are filled with shocks as described in Section [3] (see Figured]). 
The spatial diffusion coefficient are different inside and outside the bubble as described 
by Equation f[T^ . Due to the very turbulent condition inside the bubble, we consider a 
constant spatial diffusion coefficient and adopt Db = 2.08 x 10'^° cm^ s~^ (cf. estimate 
value by the fitting process in Section |5]). Outside the bubble, we take into account 
the energy (or momentum) dependence of the spat ial diffusion coeffi cient and adopt 



Dg = Do(W4GeV)0■^ Do = 6.2 x lO^^ cm^ s'^ (cf. 



Jones et al. 



200l|) 



According to the analysis in Section HI the acceleration of energetic particles in the 
bubble is facilitated by stochastic acceleration (second-order Fermi acceleration). Assuming 
that there is little or no stochastic acceleration outside the bubble, we model the momentum 
diffusion coefficient as a step function as in Equation ( IT6|) and adopt hbH^ /Db = 1-9 (i.e., 
hb = 4.4 X 10~^^ s~^ or the corresponding acceleration time scale is 7.6 Myr). 

The Galactic disk contains SNRs. We adopt the distribution suggested by 
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Stecker fc JonesI fll977l ) and modified it with a Gaussian thickness profile 



Qsnr(p, z) oc 



_P_ 



1.2 



3.22p\ 



(39) 



where p is the galactocentric radius and z is the distance perpendicular to the mid-plane. 
Here we take h = 100 pc, Rq = 8 kpc. We adopt the idea that SNRs inject energetic 
particles in the form of a power law with a high -energy cutoff at PmaxC ~ 3 x 10^5 eV. 
Therefore, together with the SNR distribution (Equation f l39|) ). the source function is 



QiP,z,p) = Qo 



-^y\xp(-^) f^V'exp (.^J^) exp (-^ 

Pma.x / \ Pmax 



(40) 



also 



As mentioned in Section El we take the SNR injection spectrum to be p = 4.35 (see 



Biermann fc Strom 



19931 ). The normalization Qo is obtained by fitting the simulation 
result to the observed spectrum and the value is 1.5 x 10^^ particles s~^ kpc~^ (GeV/c)~^. 
Integrating Q{p,z,p) over the Galaxy and momentum (from 1 GeV c~^ to Pmax) gives the 



total 
(e.g., 



uminosity of CRs 4 x 10^° erg s ^, which is consistent with the value in the literature 



Berezinskii et al. 



1990 



chap. 1, Section 4). 



Finally, the appropriate boundary conditions for the momentum coordinate are 



pdf_ 
f dp 



-4.7, /I 



P=Pup 



0, 



(41) 



P=Plow 



where the energy of the lower momentum boundary is piowC = 10^^ eV, and the upper 
momentum boundary is p^pC = 3 x 10^*^ eV. The condition at the lower momentum ensures 
that the spectral index matches that of low-energy CRs (say E < 10^^ eV). 



The spatial boundary conditions are 



dl 
dp 



0, /I 



P=PG 



0, 



(42) 



p=0 



where the radius of the Galactic disk was taken to be pc = 20 kpc. 

= 0, fU=±H = o, 



dz 



(43) 



z=0 
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where the height of the halo H = 10 kpc. 

The spectrum evaluated at Earth's position is the solid line shown in Figure O The 
model fits the data reasonably well and it is not coincident that the spectra join smoothly 
at the knee. For reference, in Figures |6] and [7] we show the spatial distribution of CRs from 
our simulation. The labels "Low" and "High" refer to the number density of particles in 
the energy range 1 x 10^^ ~ 3 x 10^^ eV and 3 x 10^^ ~ 3 x 10^^ eV, respectively. Figure [8] 
is a contour plot of the number density distribution of re-accelerated CRs {E > 3 x 10^^ 
eV) (thick lines). In the figure we also plot the distribution of seed particles from SNRs 
(thin lines). We point out that the spatial distributions of seed and re-accelerated CRs in 
the disk are quite d ifferent. In principle , CR d istribution can be derived from gamma-ray 



data. For instance, 



Breitschwerdt et al. 



(120021 ) used the gradient of gamma-ray emissivity 



in the disk derived from the EGRET data for the analysis of CR propagation in the Galaxy. 
If the diffuse gamma-ray data at E > 10^^ eV were available, the gradient test would be a 
nice tool to investigate possible proton sources in this energy range and might lend support 
to our model. 



7. DISCUSSIONS AND CONCLUDING REMARKS 

We have summarized our current understanding of the origin of CRs. It is generally 
believed that most CR power can be provided by SNRs. However the CRs with energies 
E > 10^^ eV are quite difficult to achieve in SNe due to the limited acceleration time 
and energy content in SN shocks. On the other hand, we argue that shocks in the 
Fermi Bubble produced by stellar capture events can have a much longer lifetime > 10^ 
yr and larger energy content ~ 3 x 10^^ erg, which allow them to produce CRs with 
energies E > 10^^ eV. If processes of CRs which escape from the Galaxy are taken 
into account, the predicted CR spectrum contributed by the bubble is E~'^, where 
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! + io, 



(J)/3xl030 cm2s-i)2 



3 for 10^^ eV < ^ < 10^^ eV. However, it is 



2 ' A/ 4 ' ""(m/IOS cms-i)2(/i'/10kpc)2 

very difficult to predict the exact value of the spectral index v due to the poorly measured 
value of the diffusion coefficient in the bubble. So we ffi the observed CR spectrum 
between 10^^ eV and ~ 10^^ eV and find that the diffusion coefficient is about 3 x 10^° cm^ 
s~^. By matching with the observed flux in the knee region we find that the conversion 
efficiency from shocks in the bubble to CRs is about 10%, which is quite consistent with 
numerical simulations. Other input parameters in this model, such as the capture time scale 
Tcap ~ 3 X 10"^ yr, the mean rate of energy release ~ 3 x 10^° erg s~^ per capture and 
the injected wind speed w ~ 10^ cm s~^, have b e en estimated and u sed in oth er observed 



phenomena in GC (e.g. 



Cheng et al. 



2006 



2007 



Dogiel et al. 



2009a b 



201lh . 



We put forth the idea that the Fermi Bubble acts as the re- acceleration site for the 
CRs produced by SNRs in the Galactic disk. The re-accelerated particles form the part 
of the observed spectrum that is beyond the knee (about 3 x 10^^ eV). The part within 
the knee is formed by the CRs produced by the SNRs. We demonstrated this idea by a 
numerical model. We solve the stationary transport equation in our Galaxy with two Fermi 
Bubbles (see Section [3] for our model of the bubbles). The re-acceleration process in the 
bubble is facilitated by stochastic acceleration. Our model simulated the observed spectrum 
nicely. Therefore we consider that this model provides a natural explanation of the flux, 
spectral index and matc hing at the knee of CRs in this energy range. In a related issue. 



Mertsch fc Sarkarl (120111 ) showed that the gamma ray from the Bubbles can be produced by 



stochastic acceleration of electron throughout the Bubbles. 

As described in Section [3] there are many shocks propagating in the Fermi Bubbles. 
After being re-accelerated inside the bubbles by the multiple shocks, protons (and nuclei) 



escap e the bubbles. The h 
(e.g., 



Crocker fc Aharonian 



' etime of protons by pp collision is of the order of 10^ ~ 10^° yr 



201ll ). The diffusion time for protons at these energies to escape 
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the Galaxy is of the order of 10^ yr. Thus after leaving the bubbles, protons (and nuclei as 
well) diffuse throughout the whole Galaxy (including the Earth) without any attenuation 
of energy. Not only can protons (and nuclei) be accelerated by the multiple shocks, but 
so can electrons. Nevertheless, energetic electrons lose energy efficiently. The lifetime of 
electrons can be estimated by Tg = 1/PeE (e.g., inverse Compton and synchrotron). Taking 
/3e = 3 X 10~^^ eV~^ s"^, the lifetime of electrons of energy 0.1 ~ 1 TeV is 1 ~ 0.1 Myr. 
Once they leave the acceleration site (i.e., Fermi Bubbles), they lose most of their energy 
within a short distance (less than 1 kpc). Hence the electrons are mostly confined in the 
bubbles. 



SuetaL 



(120101 ) suggested that the electron spectrum must be i?"" for E < TeV and 
the spectral index a ~ 2.4-2.8. It is clear that the energy loss processes for protons and 
electrons are very different, therefore protons and electrons in the bubble can have different 
spectra. In the CCDKI model we have assumed that the shocks can produce an injected 
electron spectrum ~ E~'^, in which the spectrum is modified by processes of energy losses 
and escape. In our subsequent work we will take into account the stochastic acceleration 
processes (multiple-shock) to see how the electron spectrum is affected. We will also further 
an alyze the spatial spectra l distribution of electrons. Preliminary results can be found 



m 



Chernyshov et al. 



(120111 ). in which the spatial distribution of gamma=ray emission is 



reproduced nicely by the multi-shock model. 
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Fig. 5. — CR spectrum at the Earth as a combination of the SNR contribution (in the 
Galactic disk) and the stochastic acceleration in the Fermi Bubbles. In addition to data 



from exper i ment 



(lAhn et al. 



s presente d in Figure HI we add ed experiments 



20081), Proton ( IGrigorov et all Il97ll ) and RUNJOB (lApanasenko et al.lboOll ). 



or lower energies: ATI C 



The black solid line is the spectrum from our numerical model. In this model, Db = 
2.08 X 10^° cm^ s~^ inside the bubbles and Dg = 6.2 x 10^*^(^0/4 GeV)^-^ cm^ s"^ outside, 
kbH^/Db = 1.9, and the injection spectrum from SNR fi = 4.35. 
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Fig. 6. — Vertical distribution of the number density of CR at two galactocentric radii, 0.4 
kpc (thick hues) and 8 kpc (thin hues), in two energy ranges, "Low" for 1 x 10^^ ~ 3 x 10^^ 
eV (sohd hues) and "High" for 3 x 10^^ ~ 3 x 10^^ eV (dotted hues). 
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Fig. 7. — Radial distribution of number density of CR at the Galactic plane (thick lines) 
and 1 kpc above the plane (thin lines), in two energy ranges, "Low" for 1 x 10^^ ~ 3 x 10^^ 
eV (sohd hues) and "High" for 3 x 10^^ ~ 3 x 10^^ eV (dotted lines). 
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Fig. 8. — Contour of the relative number density distribution of the re-accelerated CRs 
(-E > 3 X 10^^ eV) in the halo (thick lines). The seed for the re-accelerated CRs comes from 
SNRs in the Galactic plane (thin lines). 
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